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It is shown that the transverse Rayleigh Taylor like instability can be well stabilized by using 
elliptically polarized laser in the hole boring radiation pressure acceleration regime. The J x B 
effect of the laser will thermalize the local electrons and support a transverse diffusion mechanism 
of the ions, resulting in the stabilization of the short wavelength perturbations, which is quite similar 
to the ablative Rayleigh Taylor instability in the initial confinement fusion research. The proper 
range of polarization ratio is obtained from a theoretical model for the given laser intensity and 
plasma density. The stabilization mechanism is well confirmed by two dimensional Particle-in-Cell 
simulations, and the ion beam driven by the elliptically polarized laser is more concentrated and 
intense compared with that of the circularly polarized laser. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw, 52.59.-f 



Recently, ion acceleration from the interaction of ultra- 
intense laser pulse with plasmas has attracted wide at- 
tention because of its broad applications, including pro- 
ducing high energy density matter, ion- fast ignition in 
laser fusion, tumor therapy and radiographing[l|-[6|. Al- 
most all of these applications call for a high quality ion 
beam with large particle number, sharp energy spread 
and low divergence angle. Radiation pressure accelera- 
tion (RPA) is an potential scheme for generating high 
quality ion beams. According to the target thickness, 
usually there are two modes of RPA acceleration mecha- 
nism: light sail (LS) RPA for thin target 0-[i3 and hole 
boring (HB) RPA for thick target |15l423|. In particular, 
the HBRPA owns the intrinsic property for large parti- 
cle number acceleration [2^. In HBRPA, the pondero- 
motive force drives the local electrons inward, resulting 
in a shock like double layer (DL) region with large elec- 
trostatic charge separation field. The latter could trap 
and reflect the ions initially located ahead of the DL, 
compressing and accelerating them like a piston. For a 
usually circularly polarized (CP) laser driven HBRPA, 
the DL oscillations would broad the energy spread of the 
accelerated ion beams[18]. Wu et al. proposed to use el- 
liptically polarized (EP) laser to suppress the DL oscilla- 
tions, generating high quality mono-energetic ion beams 
compared with that of CP laser [25]. However the crucial 
issue, which is of fundamental influence to the RPA ac- 
celeration scheme, is the transverse Rayleigh Taylor like 
instability (RTI) [lo|, Q Hsl • The classical RTI can occur 
when a light fluid pushes or accelerates a heavy fluid, and 
this situation is quite similar to the RPA case, where the 
photons act as light fluid and plasmas as heavy fluid [14]. 
The RTI will break the target surface and terminate the 



acceleration process. However, unlike the classical RTI, 
the short wavelength ablative RTI in initial confinement 
fusion (ICF) research can be stabilized due to thermal 
smoothing of the perturbation or the transverse diffusion 
mechanism |2 7H3 1\ . 

In this paper, we propose to use EP laser for stabiliz- 
ing the transverse RTI in the HBRPA. The idea is quite 
similar to the stabilization mechanism for ablative RTI 
in the ICF research. Because of the J x B effect, the 
EP laser will thermalize those electrons located within 
the DL region, and the high plasma temperature pro- 
vides a fast transverse diffusion velocity (~ ^JT^^Z jmi) 
of the ion particles, where Tg is the plasma temperature, 
Z is the ion charge number and is the ion mass. It 
is this transverse diffusion of the ions that stabilize the 
transverse ablative RTI. During the characteristic time 
of the RTI, the diffusion range of the ions can overshoot 
the instability wavelength if a faster diffusion velocity is 
given, equally a higher plasma temperature or smaller 
polarization ratio a — a^lcLy for EP laser. However, as 
expected, if a is too small (a = for linearly polar- 
ized laser), the laser piston structure is totally destroyed. 
Thus, there should be a lower limit for the polarization 
ratio OL to sustain the HBRPA process. Based on these 
ideas, the proper range of polarization ratio is obtained 
through a theoretical model for given the laser intensity 
and plasma density. This scheme is well confirmed by 
two dimensional (2D) Particle-in-Cell (PIC) simulations. 

In the HBRPA regime as shown in Fig. 1 (a)0, IT7| . 
the ponderomotive force drives the local electrons inward, 
resulting in a shock like DL region with large electrostatic 
charge separation field. The latter could trap and reflect 
the ions, compressing and accelerating them like a piston. 
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FIG. 1. (color online) (a) Schematic structure of the piston 
maintained by the radiation pressure in the piston-rest frame, 
(b) The stabilization mechanism of the EP laser driven abla- 
tive RTI. The EP laser thermalizes those electrons, providing 
a fast transverse diffusion velocity of the ion particles, which 
stabilizes the transverse ablative RTI. 



Because of the smaller mass of electrons, the accelerated 
ion beam is accompanied with an electron beam, keeping 
almost quasi-neutral. If we think of the ion and electron 
motions as a whole, it is reasonable to neglect the electron 
inertia, and assume that the ponderomotive force acts on 
the ions directly. As the acceleration process is limited 
within the DL region, which is about tenth of laser wave- 
length, this situation can be greatly simplified sH-Q . 

Let us consider in the piston-rest frame, as shown in 
Fig. 1. The mass density of the thin layer is 
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where D[ is the DL width, n[ is the ion density distri- 
bution within the DL region and ujpi = ^^ATrZe'^ne/mi 
is the characteristic frequency of the ions. The radiation 
pressure of the EP laser is 



p,,, = 2/o(l-/3/)/(l + /3/)/c. 



(2) 



where Iq = {ay + al)ncmeC^ ay and az are 
the normalized laser amplitudes ay — eEy/rrieUjQC^ 
az = eEz/rrieUJoc^ /3f is the propagation velocity of 
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FIG. 2. (color online) The proper range of polarization ratio 
vs. laser intensity. Here, the laser wavelength is 10.0 fim and 
the hydrogen plasma density is 20nc which is 2.2 x 10^^ /cm^. 
The black line represents the upper limit derived from Eq. 
(14) and the red line is for the lower limit from Eq. (15). 
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the laser piston (3f = ^/Io/ni{mi + Zme)c^/ [1 
^/ Iq / ni{mi + Zme)c^], rrie is the electron mass, 
the ion mass, Z is the ion charge number and the term 
(1 — /3/)/(l + /3/) is the modification due to the Doppler 
effect in the piston-rest frame. The acceleration g of the 
inertial force acted on the thin layer of ions can be ex- 
pressed as 



" Prad/ ^71 



(3) 



We consider two points (xq, ^o) and (^o, ^o+^^o) on the 
thin ion layer. These two points will evolve at some time 
to the points (x, y) and {x + Syodx/dyo, y + Syody/dyo). 
The X and y components of the force equation of the 
element on the thin layer can be written as 

dpx/dt = -gamdyo + Praddyody/dyo, (4) 
dpy/dt = -Praddyodx/dyo, (5) 

where px = jfamdyodx/dt and py = jfamdyody/dt. Af- 
ter some arrangement, Eq. (4) and (5) can be rewritten 
as 



d^x/dt^ = 
d'y/dt' = 



-9hf + {dy/dyo)g/-ff, (6) 
-{dx/dyo)ghf. (7) 



The solution of Eq. (6) and (7) turns out to be[^ 

x = So exp[t(/c^/7/)^/^] cos(%o), (8) 
y = yo-So exp[t(fc^/7/)^/^] sin(%o), (9) 

where is the initial disturbed length, k = 27r/Art and 
Art is the instability wavelength. If the disturbed length 
in the x direction is 1/k at sometime r , it means that 
the adjacent sections of the thin layer begin to collide 
with each other [32]. Here the transverse RTI has already 
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developed, and we define r as the characteristic time of 
the RTL From Eq. (8), we have 



(7//fe^)^/'log(l/(5o^). 



(10) 



When transformed to the laboratory frame, the charac- 
teristic time of the RTI reads, r = 7/T , 



Ai/27^/'(m/Z)3/4ni/4(i + ^^) 
27r(a2+a2)i/2(l-/3/)i/2 



1/2 



•log(l/5ofc)(ll) 



where m = mi/rrie^ n = Ue/ric^ Uc is the critical density, 
A = A^t/Ao is the instability wavelength normalized to 
laser wavelength, and r is normalized to laser period Tq. 

The factor log(l/^o^)7 to some degree, is constant, 
which can be determined from the PIC simulations. It 
must be emphasized that log{l /Sok) is not sensitive to 
^0 and k. Once the value of log{l/Sok) is chosen, it will 
keep constant for all cases of different laser intensity and 
plasma density. According to our experience of a series 
of 2D PIC simulations, we have already defined the con- 
stant factor \og{l/Sok) — 13.8 compared with Eq. (11). 
This equation demonstrates that the shorter wavelength 
perturbations grow more faster, and that strong laser in- 
tensity and low plasma density correspond to short char- 
acteristic time of instability. 

When an EP laser falls on the target, its oscillating 
ponderomotive force will thermalize the local electrons. 

I 



Considering the electron gamma factor under the EP 
laser. 



7e = \I^^Wy cos(cjot)]^ + [dz sin(cjot)]^, (12) 



and assuming 1 + + >> — a^^ Eq. (12) can be 
transfered to, using Taylor expanding formula. 



7e = 7ec + 7eos cos(2a;ot), 



(13) 



1 + a^/2 + a^/2 is the constant gamma 
(a^ — a^)/47ec is the oscillating gamma 



where 7ec = 
term and 7e 

term. The thermalization of the local electrons is dom- 
inated by the time oscillating term, in which the tem- 
perature of the local electrons can be approximated to 
rrieC^ijeos — 1)- Thus, the transverse diffusion velocity of 
ions, equally the sound speed, is Vd = \J (7eos — l)^/m, 
where is normalized to the light speed c. 

As shown in Fig. 1 (b), in order to stabilize the trans- 
verse ablative RTI, we suppose that during the char- 
acteristic time r, the diffusion range of the ions must 
overshoot the instability wavelength |26|. To smooth the 
short wavelength perturbations with wavelength A < 1, 
we have rv^ > 1. After some arrangement, it reads 



{al-al)/4^1 



al/2- 



(14) 



where r] = [27r / log {1 / Sok)]^ = 0.21. From Eq. (14), 
the minimum temperature equally the upper limit of the 
polarization ratio a is determined for the defined laser 
intensity and plasma density. 

There must be, as expected, a lower limit for the po- 
larization ratio a to keep the piston structure intact. For 
an usually laser piston structure shown in Fig. 1 (a), the 
ponderomotive force drive the local electrons inward re- 
sulting in a shock like DL region with large electrostatic 

I 
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charge separation field, which could trap and reflect the 
ions. While the strong J x B effect of LP laser will drag 
the local electrons forward to vacuum breaking the laser 
piston structure. The lower limit of the polarization ratio 
a is derived based on the assumption that these forward- 
going electrons must be stopped within the DL region. 
Balancing the (forward) kinetic and the electrostatic and 
ponderomotive potential energies of the electrons [25], we 
can obtain the lower limit of the polarization ratio a. 



1(7/ - l)/Z + y/l + (4 - a2)Vl6 - 1 < ^/^{al + 4)l3f{l - /3/)/3(l + pfhf 



(15) 



For the given laser intensity /q = {^1 + aDncTrieC^ 
plasma density ng, ion mass m and ion charge number 
Z, we can solve out two polarization ratios a — az/ay 
from Eq. (14) and (15). These two values are the upper 
and lower limit of suitable polarization ratios. For the 
hydrogen plasma with density ng = 20nc, the proper 



I 

range of polarization ratio a = az/ay can be obtained 
according to Eq. (14) and (15), which is shown in Fig. 2. 

2D PIC (KLAP code [34(33) simulations are run to fur- 
ther confirm this scheme. The size of the simulation 
box is Lz X Ly = 20Ao(2;) x 40Ao(^) with Aq represent- 
ing the laser wavelength. The simulation box is divided 
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FIG. 3. (color online) (a) and (b) The distribution of proton density driven by CP and EP laser pulse at t = 25To, t = 35To, 
t — 40To and t — 45To respectively. Here, the laser wavelength is 10.0 /xm, the hydrogen plasma density is 20nc which is 
2.2 X 10^^ /cm^, the laser intensity is 6.85 x 10^^ W/cm^ and the polarization ratio a — 0.70 for the EP laser pulse. 



into uniform grid of 2000 (z) x 4000(?/). The CP laser 
pulse enter into the simulation box from the left bound- 
ary. The bulk target consists of two species: electrons 
and photons, which are initially located in the region 
l.OAo <z < 2O.OA0 and -20.0Ao <y < 20.0Ao with den- 
sity ng = 20nc, where Uc = cjQe^me/47r = 1.1 x 10"*^^ /cm^ 
is the critical density for 10.0 /im laser pulse [23|, [ssj . 
We use 160 particles per cell to run the simulations. 
The initial plasma temperature is set to be 10.0 KeV. 
To exclude the effect like electron rebound on the target 
boundary, we applied particle absorbing boundary con- 
ditions, where the absorbed particles are compensated 
by the incoming particles with the background tempera- 
ture from the boundary to ensure charge neutrality. The 
normalized amplitude of the CP laser electric field is 
ay = 50.00 and = 50.00, corresponding to the laser 
intensity 6.85 x 10^^ W/cm^. To exclude the effect of 
target surface curvature caused by the shape of the laser 
pulse, we applied transverse super-Gaussian laser pulse. 
The laser pulse has a temporally Gaussian rising profile 
exp[— 1/25)^] followed by a constant intensity, where t is 
normalized to laser period Tq. In contrast, the EP laser 
pulse with the same space and temporal profile is also 
run. According to Fig. 2, for laser intensity at 6.85 x 10^^ 
W/cm^, we choose the polarization ratio a = 0.7, which 
equally corresponding to ay = 57.93 and = 40.55. 

Our main results are clearly shown in Fig. 3. When 
driven by CP laser, the ponderomotive force has no 
second-harmonic oscillating component, and there is no 
thermalization effect at all. As shown in Fig. 3 (a) 
the transverse RTI has rapidly developed at t = 40To, 
and the perturbations turn out to be even more fierce 
at t = 45To, with the perturbation wavelength about 



1 ~ 2Ao. These short wavelength perturbations can 
severely disturb the hole boring process, terminating the 
HBRPA acceleration at early time. In contrast, when 
driven by EP laser, because of the J x B effect, the EP 
laser will thermalize those electrons located within the 
DL region, and the high plasma temperature provides 
a fast transverse diffusion velocity (~ ^JT^ZJmi) of the 
ion particles. In this case, the temperature of the lo- 
cal thermalized electrons can be as high as 3 MeV, thus 
the diffusion velocity can reach as fast as 'U^i = 0.05. It 
is this transverse diffusion of the ions that stabilize the 
transverse ablative RTI, which is clearly shown in Fig. 3 
(b). During the characteristic time of the RTI, r ~ 40To, 
the diffusion range of the ions can be as far as 2.0Ao. 
The range of the transverse diffusion ions can readily 
overshoot the perturbations of short wavelength about 
1 ~ 2Ao. In addition, as expected, the laser piston struc- 
ture is kept intact. The HBRPA process sustains for a 
rather long time compared with that of CP laser. 

The energy spectrum of the proton beam driven by 
the EP laser is much better than that of the CP laser, 
because of the relatively smooth laser piston surface or 
shock like electrostatic charge separation field, which is 
clearly shown in Fig. 4. For the EP laser, because of 
the stabilization of the transverse RTI, the number of 
the accelerated proton particles is always accumulating 
resulting from the stable HBRPA process. While for the 
CP laser, the accumulating process is terminated at t = 
40To, and the chaotic heating is performed instead. 

Fig. 5 shows the angular distribution of the acceler- 
ated proton beams by CP and EP lasers respectively: 
(a) is for CP laser and (b) for EP. Compared with (a) 
and (b), we conclude that the transverse diffusion mech- 
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FIG. 4. (color online) (a) and (b) The energy spectrum of 
the proton beams driven by CP and EP lasers at t = 25To, 
t — 35To, t — 40To and t — 45To respectively. Here, the laser 
wavelength is 10.0 /xm, the hydrogen plasma density is 20nc 
which is 2.2 x 10^^ /cm^, the laser intensity is 6.85 x 10^^ 
W/cm^ and the polarization ratio a — 0.70 for the EP laser 
pulse. 




FIG. 5. (color online) (a) and (b) The angular distributions 
of the proton beams driven by CP and EP lasers respectively 
at t 45To. 



tion of the short wavelength perturbations, which is quite 
similar to the ablative Rayleigh Taylor instability in the 
initial confinement fusion research. The proper range of 
polarization ratio is obtained from a theoretical model for 
the given laser intensity and plasma density. The stabi- 
lization mechanism is well confirmed by two dimensional 
PIC simulations, the ion beam driven by elliptically po- 
larized laser is more concentrated and intense compared 
with that of circularly polarized laser. 
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anism driven by EP laser does not have a obvious affec- 
tion on the angular distribution of the proton beam. The 
transverse diffusion velocity of protons in this situation 
is around Vd = 0.05, while the longitudinal velocity is 
double the piston forward velocity, which can be as high 
as 'U/ = 0.4. It is reasonable to say that the small di- 
vergence angel can still be held under the mechanism we 
proposed. It demonstrates that the divergence angels of 
proton beams driven by both CP and EP lasers can be 
maintained within 10 degree. However, the proton beam 
driven by EP is more intense and concentrated than that 
of CP. 

It should be emphasized that different polarization ra- 
tios with 0.65 < a < 0.75 are also run. The stabilization 
mechanism is not quite sensitive to the polarization ratio. 
The scheme is confirmed to be rather robust. 

In summary, we propose to use EP laser to stabilize 
the transverse RTI in the hole boring radiation pressure 
acceleration regime. The J x B effect of the laser will 
thermalize the local electrons and support a transverse 
diffusion mechanism of the ions, resulting in the stabiliza- 
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